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Abstract 

We investigate the propagation modes of gauge fields in an infinite RandaU-Sundrum scenario. 
In this model a sine-Gordon soliton represents our thick four-dimensional braneworld while an 
exponentially coupled scalar acts for the dilaton field. For the gauge-field motion we find a differ- 
ential equation which can be transformed into a confiuent Heun equation. By means of another 
change of variables we obtain a related Schrodinger equation with a family of symmetric rational 
(7 — ioz'^)/{l — z^)^ potential functions. We discuss both results and present the infinite spectrum 
of analytical solutions for the gauge field. Finally, we assess the existence and the relative weights 
of Kaluza-Klein modes in the present setup. 
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I. INTRODUCTION 



One of the main purposes of superstring theory is the inclusion of all the relevant fields 
of Nature together in one single Lagrangian. Field theoretic scenarios inspired in such a 
theory put in contact gauge and matter fields with metric degrees of freedom, altogether 



defined in some extra-dimensional space Extra-dimensions, combined with the influence 
of the gravitational field, modify nontrivially all sectors so that gauge forces must be proven 
to remain the same in the usual four-dimensional (4D) subspace or predict new physics in 
some consistent way. Localization in the gauge sector is expected to hold and the effective 
4D electromagnetic force must be mediated by massless photons as usual. On the other 
hand, higher- dimensional spaces help solving fundamental problems such as the hierarchy 
gap between the Planck and gauge-coupling scales in the Standard model 

Also inspired in String theory is the use of branes to represent our Universe. In String 
theory, gauge modes are deposited on D-branes from open strings ending on them, so we 
expect gauge fields in field-theoretic models to have finite localized modes on stringy topo- 
logical defects of lower dimensionality. Actually, to obtain finite eigenstates in 4D it has 
been shown that to fulfill this task we need not just gravity but also a dilaton js, 4|, a field 
already predicted in String theory. 

In the present paper, we describe gauge fields in a warped five-dimensional bulk with 
a dilaton and a brane defect that mimics the ordinary world. Both brane and dilaton 
configurations are geometrically consistent solutions of a two scalar world action in a curved 
5D space-time where the field potential is of sine-Gordon type. 

We show a relationship between the fundamental parameters of the 5D theory which is 
crucial to determine the dynamics of the fields both in the bulk and ordinary space. Indeed, 
for different choices of a parameter defined by the quotient of some power of the sine-Gordon 
frequency-amplitude and the 5D Planck mass, the equations of motion of the gauge field can 
be completely different. Notably, in Ref. ^| we have been able to find the whole spectrum 
of a theory involving both Maxwell and Kalb-Ramond fields for a particular value of this 
parameter. As we will see, there exists a minimal value for the dilaton coupling constant 
above which the finiteness of the action is assured and it is directly related to the localization 
of gauge fields. 

In what follows we analytically obtain the propagation modes (massless and massive) of 



a gauge theory in a background of the sine-Gordon type that results in new equations of 
motion. We show that the dynamics of the quantum mechanical system associated with 
the problem is given by a simple (rational) potential function and that the solutions to the 
Schrodinger equation are of the Mathieu type (with a power-law factor). In a more general 
case we obtain the exact spectrum given by the set of confluent Heun functions and show 
that Kaluza-Klein states are strongly suppressed in ordinary space. 

The paper is organized as follows. In the next Section, we present the geometrical back- 
ground. In Section [Till we introduce the action for the 5D gauge field coupled to warped 
gravity and a dilaton background and derive the 5D equations of motion. In Section |IV] we 
obtain the quantum analog problem showing explicitly the quantum-mechanical Schrodinger 
potential. The eigenvalue spectrum is computed and graphically shown. Next, in Sections 
IVl and I VI I we discuss the general problem and draw our conclusions. Other recent results 
about thick braneworlds can be found in e.g. jof. 



II. GEOMETRICAL BACKGROUND 



Our framework is a five-dimensional space-time embedding a four- dimensional membrane 
also called thick brane. The (space-like) extra-dimension is assumed infinite and the brane 
will be dynamically obtained as a solution to the Einstein equations for gravity coupled to 
a pair of scalar fields. One of these scalars represents a domain wall defect (the thick brane) 
while the other is the dilaton. The dilaton, together with the warping of the fifth dimension, 
happens to be crucial in the gauge theory that will be developed and makes more clear 
the stringy origin of the theory. Since gauge field theory is conformal [7] all the information 
coming from the warping of the 4D metric is automatically lost. As a consequence the photon 
is non-normalizable in the four-dimensional space unless the gauge coupling is dynamically 
modified. Indeed, the exponential coupling of both the dilaton and the 5D warping to the 
gauge field conveniently modifies the scaling properties and the zero-mode becomes localized 
3|. 

The five- dimensional world action which determines the background is 

Sb = J d^x dy^/- det Gmn [2M'R - ^{d^Y - ^{dUf - Vi<^M, (1) 
where M is the Planck mass in 5D, and R is the Ricci scalar. The solution for $ represents the 



world membrane and the corresponding field solution for 11 will be the dilaton configuration 
consistent with the metric and the kink. As usual we adopt Latin capitals on the bulk and 
Greek lower case letters on 4D. 

We next adopt the following ansatz for the metric 

ds^ = e^^^y^r]^,dx^'dx' + e^^^y^y^ (2) 

where A and S depend just on the fifth coordinate, and diag(r7) = (—1,1,1,1). The 
equations of motion for eq.([T]) are 

^($0' + ^(n')' - e2^(^)V(<l>, n) = 24M3(A')', (3) 

+ ^(n')^ + e2^(^)V($, n) = -12M^A" - 24M=^(A')^ + 12M3A'S', 



and 



$" + (4A' - = e^^ ^, 

a$ 

n" + (4A' - s')n' = e'^ ^. (4) 

aU 

where the prime means derivative with respect to y. 

By means of a supergravity motivated functional >V($) defined by 

--'^ 

js] the system of differential equations can be more easily handled. This method is also 
applicable to non-supersymmetric domain walls jo], [lo| as the present one. 

First, we consider the action in the absence of gravity (and no dilaton) in order to obtain 
an expression for $. Then, we put this solution into the equations of motion ([3]) and (|4]). 

The standard sine-Gordon Lagrangian reads 

LsG = -\d^^ - V{^) (6) 

with 

y($) = ^ (1 -cos(6<l>)). 

The free parameter h signals bulk-symmetries 5$ — )■ 2n7r/6 {n G Z) among the vacua of 
this theory. Solutions interpolating vacua are possible and, assuming they depend only on 
y, one-solitons read 
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$(?/) = - arctan e^. (7) 



These functions kink on our 4D-world slice, namely at ?/ ~ 0. 

In a gravitational background of the form ([2]), now including also the dilaton, the equa- 
tions of motion and (jl]) are still compatible with solutions ([7]) provided we find the 
appropriate potential functional V for the general action ([T]), viz. 

V($, n) = exp {U/VUAP) l^li^r - • (8) 

Taking into account eq.(l5]), the superpotential functional results 

W($) = -lcos(^$) (9) 

and then 

V($, n) = -e(n/v^) ^ $ + -A-^ eos^ ^ $) . (10) 

If we now conveniently write the Hamiltonian a la Bogomol'nyi, we can detect the following 
relations among the warping functions, the dilaton and the superpotential 



U = -V3APA, S = A/4, A' = -W/12M^. (11) 
Finally, totally solving the equations of motion, the dilaton field is given by 

and 

A = -^lncoshy, S = -Y^l^^°^hy. (13) 
The relation between 11 and $ allows also writing V as 

V(*, ^ 4 (.„ ^*)./««'^' (l + - 1) ™^ ^ «) . (14) 

which fully shows its dependence on b and M (see Fig. [1]). 

As it happens with dilaton configurations related to D-brane solutions, functions such as 
([7]) and (1121) are singular when \y\ — t- oo. However, since the metric vanishes exponentially 
and both dilaton and warp factors operate under an exponential coupling, the model is kept 
free of divergences. 

The warping functions amount to a shift in the effective four-dimensional Planck scale, 
which remains finite with the following definition 

/oo 
dy e4^(^)+^(^). (15) 
-oo 
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Figure 1. Family of background potential functionals V($) (eq ]14p for different values of a = 
even a = 2,4,6,8,10 (solid line), odd a = 1,3,5,7,9 (dashed line). For clarity we 
adopted thinner lines for bigger values of a. 

Using the consistency relation (fTTl) just found, the action reads 



S 



B 



(4) 



(16) 



where S'(4) is the remaining of the action integrated in 4D. According to the solution 
A(y) = 2alnsechy (c.f. eg. (Ill p - eq.(IT3l)) the 5D factor results finite provided c = (17 -|- 
2\VMP)/4: > 0, namely A > = Aq. 

Studying the fluctuations of the metric about the above background configuration, it 
is possible to see that this model supports a massless zero-mode of the gravitational field 
localized on the membrane even in the presence of the dilaton. In order to prove the stability 
of the background solution, we would have to show that there are no negative mass solutions 
to the equations of motion of a perturbation /i^i, of the metric. Actually, a gravitational 
Kaluza-Klein spectrum appears, starting from zero and presenting no gap. This can be 
easily seen after an appropriate change of variables and decomposition of the gravitational 
field, and a subsequent supersymmetric type expression of the Schrodinger type operator 
resulting from the equation of motion (see js, 11 1 for details). The issue of the coupling of 
these massive modes to the brane has been analyzed in detail in jl2| . 



III. GAUGE FIELD ACTION IN A WARPED SPACE WITH DILATON 



Let us consider the following 5D action where a five- dimensional electromagnetic field 
An is coupled to the dilaton [l3| in a warped space-time 



Sg = J dyd^x^Z-detGAB e"^" j-^FM^F^^^j 



(17) 



where Fmn = diM^N]- 

Assuming that the gauge field energy density should not strongly modify the geometrical 
background, we can study the behavior of the propagating modes in the background of the 
topological configuration studied in the last Section. In general, most of the attempts to 
stabilize 5D brane worlds by means of a scalar field in the bulk do not take into account the 
back-reaction of the scalar field on the background metric js], 11, Q, 14 1 and those in order 
to conipute the scalar back-reaction on the metric were unsuccessful except in a few special 
cases jol. Il5|. 



The factor exp(S — An/2), present in the integrand, will lead to a change in the integration 
measure which is crucial to conserve the effect of the warping function on the 4D gauge field. 
As a consequence, zero-modes (namely photons) are normalizable in the 4D effective theory. 
To see it in detail, we need to solve the 5D equations of motion for Am 
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where diag Gmn = (e^^?7/^i/, e^^). For this, we adopt the following gauge choice A^ 
0, dfj^A^^ = and separate the fifth from the ordinary coordinates as follows 



(18) 



Now, from eq. llTSil we just get 



A>'{x,y) = a^'ixHy). 



1 



(19) 



[n + ^d,{fd'u)] a^ = 0. 



(20) 



Note that the warped metric and the dilaton deform the solutions of this differential equation 
by means of the factor f(y) = e4A+s-An/2 multiplying u{y) and u'{y). A full Kaluza-Klein 
spectrum results from the solution of the general case 



(21) 
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where is an arbitrary constant representing the 4D squared boson mass of the vector 
gauge field. It means that = a^(0)e*'P^' with p"^ = —m?. 

By expanding eq. (12U]) . we obtain the most general y-dependent equation of motion for 
the modified sine-Gordon potential f ll4l) derived from action f[T71) in a form which exhibits 
its dependence on a = 1/6M^6^, and c 

u"{y) + a(l — 2c) tanhy u\y) + sech" y u{y) = 0, (22) 

where y G (— oo, oo) as already stated. Looking back at the definition of the auxiliary 
constants we get the explicit dependence of the solutions on the original parameters b, A 
and M. 



Below, we will discuss the possible values of as resulting from an eingenvalue problem 
related to the equation of motion (!22l) . Indeed, there exists a Schrodinger like equation 
equivalent to eq. fl2^ with a potential function which concentrates all the richness implicit in 
the complicated equation ( l22i) . Note that the particular solution u{y) =constant represents 
the = photon state of the 5D theory. Since this solution satisfies eq. (l22!) for any value 
of a and c, any member of the family of problems has a guaranteed localized zero-mode. See 
[sl for details. 

Localization of gauge-field modes in the ordinary space can be established by verifying 
that the corresponding 5D action is finite. From eq. (fT9ll one has F'^'^ = f'^'^u{y), where 
j!fj.u _ QfiaQup j^^^ gQ ^YioX, for a gauge mode A^^- the relevant part of eq. lITTl) reads 

SA^t] = jdy «2(^)e4Afe)+^(^)-An(.)/2 j (23) 

Using the field solutions found in eq.(I7]) and the equations thereafter, the fifth dimension 
factor will remain finite for each mode uiy) growing below e'^'^ at infinity. Thus, any finite 
solution is a physically acceptable Kaluza-Klein state (as we have seen above, to have a 
finite 5D Planck mass and background action Sb we already need c > 0, i.e. A > Aq). 
It is known that by means of a transformation 

fly 

u{y) = e-'^'^/^Uiz), ^ = e-^^ (24) 
we can turn eq. (l22|) into a Schrodinger-like equation in the variable z (see e.g. 

0,3). 

In 

general, the existence of an analog Schrodinger equation is useful to give us a feeling of the 
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physical profile of the solutions of the original problem, as e.g. parity and eigenvalues. With 
a = c — 1/4 and /3 = — 1/4 we can eliminate the first derivative term in U and have a pure 
mass term as usual. The resulting differential equation reads precisely 



U{z) = m'U{z), 



where ^aiz) = e-^/^j^f A" - -fA'^) and 7 



In a few cases the last expression 



can be inverted after exact integration in order that an analytical expression for the analog 
non-relativistic potential comes about. In Ref.jsl we have solved the a = 2 case and found 

QJ2(z) = -2a [1 - {2a - 1) tan^ z] . 



In this paper, for a = 4 we find 



(7 - uz^) 



(1-^2)2 ' 

where 7 and u are constants and we shall analyze it in what follows. 



IV. THE QUANTUM ANALOG 



In the present case we can turn eq.f l22p into a Sturm-Liouville problem by means of 

z = tanhy, (25) 

u{y) = (coshy)^'-^ U{z) (26) 

(see eq. (l2^ ). Now, we have a related Schrodinger equation defined in the z variable 

" d^ 



dz^ 



with 



+ 2J4 u 



4c 



U{z) =m'U{z), 



[l + 2(l-2c)^2] . 



(27) 



(28) 



(1-^2)2 

We can see that the potential function diverges at z = ±1 and so the boundary conditions 
of this analog problem are {U{z = ±1) = 0, U'{z = ±1) finite} which must be in order to 
match finite u{y) solutions to eq. ( l22l) at ?/ — ?■ ±00. After solving the quantum analog we 
have to transform back variables and functions to check the finiteness and continuity of the 
original solution u{y) in order to be physically acceptable. 



We now better introduce the variable 6 by means of 



z = cos9 (29) 

which results in equation 

U"{9) - cot{9) U'{e) - ^^—^ \l + 2(1 - 2c) cos^ ^1 U{e) = -m^ sin^O U{9) (30) 

sin 6 

for the analog wave function U{6) with U{6 = vr, 0) = 0. 

According to the arguments of localization seen in the previous section, physically ac- 
ceptable solutions require c > so we shall be restricted to that region. 



1. The c=l/4 case 

Equation (l30l) gets strongly simplified for the value c = 1/4. In this case we obtain 

1 d / 1 



sin 9 d9 \ sin 9 



( -^U'(9)] + sm^9 U(9) = (31) 
\sin6' / 



with solutions 



f/(i) (9) = Uo sin(m cos{9)) (32) 
f/(2)(^) = f/ocos(mcos(^)), (33) 

which in terms of the original variable and function read 

'^^^\y) = Mo sin(mtanhy) (34) 
u^'^Hu) = Mo cos(mtanhj/), (35) 

(see Figs. U\ E])- The zero-mode, m = 0, is then related to u{y) = uq as already mentioned. 

Since 5J(z) is an even function (in this case trivial), solutions must have definite parity. 
Besides, the potential divergence at 2; = ±1 implies that the corresponding solutions are 
expected to be zero there. Thus, antisymmetric solutions correspond to the eigenvalues of the 
Schrodinger equation (1271) m = mr while symmetric solutions correspond to m = (2r2+l)7r/2, 
with n G (or simply m = {n + 1)tt /2 with n even for symmetric solutions and odd for the 
antisymmetric ones). 
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Figure 2. Plot of sin(mtanh(y)), for c = 1/4, m = tt (black line), m = 2tt (long-dashed blue line), 
and m = Stt (dashed red line)). 




Figure 3. Plot of cos(mtanh(y)), for c = 1/4, m = 7r/2 ( black line), m = 37r/2 (long-dashed blue 
line), and m = 57r/2 (dashed red line)). 

2. Other analytical solutions 



If we perform the transformation 



U{9) = sin^^ M{9), 



(36) 



in place of eq. (l30|) we get the following problem for ^A{9) 

M"{e) + {2k - 1) cot{e)M'{e)+ 



1 — Ac 

k{k - 2) cot\e) -K + —^[1 + 2(1 - 2c) cos^e] 



sm'9 



M{9) 



m'sm^OMie). (37) 



Now, we can choose a convenient power for the last transformation, k = 1/2, in order to 
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turn this into 

M"ie) - 4 



M{e) = -m'sin'OMie) (38) 



which, for 4c — 4c + y| = 0, is known as the Mathieu differential equation 



M"{9) + {Ac-3/2 + m'^sm^e) M{e) = 0, 



(39) 



with c = 3/8orc = 5/8. 
3. The case c = 5/8 



In this case, Eq. (JSHD results in 

(777 ^ 777 ^ \ 
l + - ^cos(2^)jM(^^) = 

whose analytic solutions are the general Mathieu functions 



(40) 



Table I. List of the first 21 values of rUg (symmetric solutions) and nia (antisymmetric ones) for 
c = 5/8 





•ma 


0.0000000 




4.0649860 


2.3807959 


7.2962115 


5.6914019 


10.4780880 


8.8902613 


13.6431458 


12.0619596 


16.8002756 


15.2224185 


19.95298930 


18.3770538 


23.1029720 


21.5282531 


26.2511409 


24.67724335 


29.3980716 


27.8247238 


32.5440439 


30.97112352 
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Figure 4. Symmetric and antisymmetric eigenfunctions U in z space for c = 5/8;m = (black line, 
symmetric), m = 2.380795874 (dashed blue line, anti-symmetric), m = 4.0649860 (long-dashed red 
line, symmetric), m = 5.6914019 (dotted black line, antisymmetric), and m = 7.2962115 (solid 
green line, symmetric). 

M^'\9) = Mc(l + ^,^,e^ (41) 
M(^)(0) = Ms(^l + ^,^,^), (42) 

which, in terms of the z variable, result in the analog wave-functions U{z) 

[/(I) (z) = (1 - ^2)1/4mc (^1 + arccos(z)^ (43) 

(22 \ 
l + ^,^,arccos(;2)J (44) 



(see [l6| for details about Mathieu functions). 

As mentioned above, the boundary conditions of the present problem are {U{z = ±1) = 
0, U'{z = ±1) finite}, related to finite u{y) solutions to the original equation, recalling 
that y E (—00,00). The first set of solutions, U^^\z), is not physically interesting be- 
cause the derivatives of these functions are divergent at the boundary. The reason is that 
Mc(arccos(z)) cannot be zero at 2; = 1 for any value of m. The second set, on the other 
hand, has physically acceptable solutions for a discrete set of values of m, the (twenty) first 
of which we show in Table [H These solutions are symmetric or antisymmetric, as expected 
(see Fig. H]). Note the presence of a zero mode. 
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Figure 5. Symmetric solutions u^'^\y) (Eq. ([16|) ) for c = 5/8; m = (dash-dotted black line), 
m = 4.0649860 (long-dashed blue line), m = 7.2962115 ( dashed red line), and m = 10.4780880 
(solid black line). 
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Figure 6. Antisymmetric solutions u^'^\y) for c = 5/8; m = 2.3807959 (solid black line), m = 
5.6914019 (long-dashed blue line), m = 8.8902613 (dashed red line), and m = 12.0619596 (solid 
black line). 



In terms of we have 



u^^\y) = coshy Mc ( 1 + — , — , arccos(tanh?/) ) (45) 
u^'^\y) = cosh?/ Ms ( 1 + — , arccos(tanh?/) ) , (46) 



where the set u^^\y) diverges when y — > ±oo, as due from the comments above, so we just 
keep the solutions u^'^\y) (see Figs. [5]and|6]). 
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4- The case c = 3/8 
Now, Eq. ([MD reads 



with solutions given by 



M"{e) +(^-'^ cos(2e)^ M{e) = 



M^'\d) = Ms[—,-^,d 



1 1 
m m 



2 ' 4 ' " 



corresponding to 



z 



:i-^2)i/^Mc(^^,^,arccos(^) 



U^^\z) = (1 - z'f'^Ms ^,arccos(z 
in the z space with the boundary conditions aheady seen. 



Table II. List of the first 20 values of rus and nia for 



m. rrir, 



1.14718042 

4.30206964 

7.44879288 

10.59305044 

13.73629872 

16.87903028 

20.02145951 

23.16369547 

26.30579984 

29.44781028 



2.72632477 
5.87592007 
9.02110025 
12.16475984 
15.30771212 
18.45027434 
21.59259704 
24.65270800 
27.87681505 
31.01878823 
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Figure 7. Symmetric solutions U^'^^z) for c = 3/8; m = 1.14718042 (solid black line), m = 
4.30206964 (long-dashed blue line), m = 7.44879288 (dashed red line), and m = 13.73629872 (solid 
green line). 




Figure 8. Antisymmetric solutions U^'^\z) for c = 3/8; m = 2.72632477 (solid black line), m = 
5.87592007 (long-dashed blue line), m = 9.02110025 (dashed red line), and m = 15.30771212 (solid 
green line). 

As we discussed in the previous (c = 5/8) case, only the second set of solutions is 
physically relevant and just for a discrete (infinite) sequence of m eigenvalues. For such 
values solutions have definite parity according to 5J(-2), Eq. fj28|) . (see Figs. ([7]) and ([8])). 

Note that the solutions to the Schrodinger equation (the analytical expressions fl5Up and 
f l3T]) ) are not compatible with a zero- mode for the z-boundary conditions given above. Ac- 
tually, as a general result, for any value excluded from the sequence starting in Table [TTl 
(z) -eq. (!50|) - has divergent derivatives at z = ±1 and U'^'^^z) -eq. (15T]) - is not even 
symmetric. For this reason the zero-mass solutions of eq. ( H7|) . M{6,m = 0) G {cons.,^}, 
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Figure 9. Mathieu function at the boundary (6 = vr) as a function of the mass for c = 3/8, 
Ms (^^,^,6 = Try At the other boundary, Ms (^^,^,0 = = 0. 

do not correspond to valid solutions of the Schrodinger problem. In FigOwe can see all the 
first mass values of the sequence which nullify the Mathieu functions Ms at the boundary. 
These values, also listed in Table UTt guarantee finite derivatives of U^'^^z). The absence of 
the zero mode in this list indicates a limitation of the Schrodinger analogue approach. We 
will come again to this point in the next Section. In terms of y we have 

— , — , arccos(tanh?/) j (52) 

— , — , arccos(tanh2/) j , (53) 

of which M*^^^ represents the only non-divergent set of solutions, as we illustrate in Figs. ( ITOl) 
and (fTT]) for the first quantum values of m. 



V. THE CONFLUENT HEUN EQUATION 

We now investigate our original problem by relaxing the quantum analog condition. In 
order to obtain the general solution of Eq. fl22|) we perform the following change of variable 

X = tanhy. (54) 

It maps the y space to x G (—1, 1) and we will eventually transform it back in order to 
come into the original space and variables. Note that Eq. fl22]) is symmetric under a parity 
transformation and thus the differential equation admits even as well as odd parity solutions, 
as it should. Now, Eq. ( 122|) becomes 
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u"{x) + (2 - 5)-^ — -u'{x) + m^{l - x^)'^-\{x) = 0. (55) 

which is an homogeneous second-order hnear differential equation with polynomial coeffi- 
cients provided a is even. Here c = a(l — 2c) so that c = {—oo, a). Now, Eq. (j55l) looks more 
familiar if we change into z 

/// ^ l-c/2\ , m^u(z) ^ 
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Figure 10. Symmetric solutions u^'^\y), eq. (j53p . for c = 3/8; m = 1.14718042 (solid black line), 
m = 4.30206964 (long-dashed blue line), m = 7.44879288 (dashed red line), and m = 13.73629872 
(solid black line). 
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Figure 11. Antisymmetric solutions u^'^\y), eq. ()53p . for c = 3/8; m = 2.72632477 (solid black line), 
m = 5.87592007 (long-dashed blue line), m = 9.02110025 (dashed red line), and m = 15.30771212 
(solid black line). 
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for the case under study. Formally, this equation has two regular singular (Fuchsian) points 
at z = 0, 1, and an irregular one ai z = oo. This is known as a Confluent Heun equation 
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19|. 



We can compare Eq. fISBl) with the canonical non-symmetrical general form of the con- 
fluent Heun equation as given in 

Hc"{z) - 



a 



z-1 



+ 



[5 + f (/3 + 7 + 2)];^ + r/ + f + 1(7 - «) (/3 + 1) 



Hc{z) = 0, 



(57) 



z{z-l) 

whose solutions around z = are denoted by 

H^^^ = Hc{a,P,-f,6,r];z) (58) 

H^^) = z-f^ Hc{a,-/3,-f,6,7];z). (59) 

In general, there are two linearly independent local series solutions around each singular 
point. In the region of interest, 2; < 1, we look for a regular local solution around z = 
which is defined by the Heun series as 

00 

Hc{z) = J2dnz''. (60) 

n=0 

Here the constants dn (with d^i = and do = 1) are determined by the three-term recurrence 
relation [2l| 

Andn = Bndn-l + Cndn-2, 



where 



Br, 



P 1 
1 + - ^ 1 

. -a + /3 + 7-l +(a-/3-7)/2- a/3/2 + /37/2 



n 



^ 1 + 



-c/2-3/2 g/2 + l/2-mV4 



n 



^ I a;(/3 + 7) , m 

C„ = — (5 + \. + a n - 1) 1 ^ 



Am? 



(61) 
(62) 

(63) 
(64) 



By comparing Eqs. ( 156|) and ( |57I) . it is easy to identify a = 0, /3 = —1/2, 7 = — c/2, 
5 = m^/4, and r7 = c/8 + l/4 — m^/4. Then the solutions of Eq. f l22l) are given by 

(65) 



u^^Hy) = i/c 0, — , — , — , - + -; tanh^y 

\y) \^ ' 2 2 4 4 8 4 

('21/ \ / 1 c 1 c 2 

^(y) = tanh?//7c 10, 2'~2'X' 4 ^ 8 ~ X' * ^ 



(66) 
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for arbitrary values of c (or c), namely of the dilaton coupling constant. The conditions 
these Heun u{y) solutions must obey to be acceptable are the original ones, i.e. finiteness 
and continuity in the whole space. 

Table III. List of first values of and rua for c = —30. 





8.69355330 

13.16860126 

17.10643340 

20.80307154 

24.36269818 



rria 



5.90953031 

11.03547247 

15.17836781 

18.97642242 

22.59620420 




(a) 

Figure 12. Symmetric solutions of the Heun equation (c = —30) for m = (black dash-doted), 
m = 8.69355330 (blue dashed), m = 13.16860126 (red dashed), m = 17.10643340 (black solid), 
m = 20.80307154 (brown solid), m = 24.36269818 (orange solid). Curves required masses with 
thirty significant digits of which only the first are shown, (a) Solutions eq. (|65p near the origin 

A noteworthy point in the present approach is that now, depending on c, the mass values 
can be quantized, as we saw in Sect. IIV[ or not, as we will explain in what follows. 

After a lengthy numerical exam, we found clear evidence that for c < (namely A > 
Ai = 15/17Ao) all the mass spectra are discrete. For c G (0,4), Aq < A < Ai, on the other 
hand, the corresponding spectra start with a zero mode and grow continuously. This sharp 
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(a) 

Figure 13. Antisymmetric solutions of the Heun equation (c = —30) for m = 5.90953031 (black 
solid), m = 11.03547247 (blue dashed), m = 15.17836781 (red dashed), m = 18.97642242 (brown 
solid), (a) Solutions eq. ([S5|) near the origin. 

contrast may be traced back to eq. (1221) where the second term of the differential equation 
flips precisely with the sign of c. Note that for any well-behaved solution u{y), the third 
term of eq. (1221) can be disregarded at infinity. The remainder differential equation can be 
easily solved showing that, for c > 0, solutions are always convergent to zero and for c < 
they diverge at the boundary. Guided by this result, we performed a numerical survey in 
each region arriving at the conclusion above: for c > there exist physical solutions for 
arbitrary m while, otherwise, only a discrete sequence of masses allow for finite solutions at 
the border. 

It should be mentioned that for small values of c the solutions stabilize quickly. On the 
contrary, for c < —10 the numerical calculation is more difficult and more digits are needed 
in the mass precision to stabilize solutions at large values of y. For example, for c = —30, 
which corresponds to A = 0, more than thirty significant digits were necessary in the mass 
spectrum to find the solutions as shown in Figs. [I2] and [131 In Table IIIII we listed the first 
values of the mass up to the eighth decimal place. 

As awaited, for the cases studied in the previous Section we find again the same results. 
Note however that the zero-mode in the c = 1 (c = 3/8) case now appears explicitly. 
This was expected since there exists an analytical m = solution to Eq. ( 122|) . namely 
uo{y) = ci arctan(e^) + 62, which must be present in a full approach. Furthermore, for c = 1 
the Heun solution indicates that Table II would not only start from zero but would also 
be continuously filled in as mentioned above. In Figs. [T3] and [15] we can see finite analytic 
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Heun solutions, given by eq. ( 165|) and eq. ( 166|) . for some arbitrary values of m besides the 
quantum-mechanical analog ones. Another way to see it is by means Fig. [16] and Fig. [T7] 
where m has been fixed arbitrarily to one of the eigenvalues of c = 1 and c is then varied. 
In the c = — 1 (c = 5/8) case the spectrum still obeys quantized values as given in Table I. 
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Figure 14. Symmetric solutions Eq. (j65p for c = 1; m = (dash-dotted black line), m = 1.14718042 
(solid black line), m = 4 (dot green line), m = 4.30206964 (solid blue line), m = 7 (dashed black 
line), and m = 7.44879288 (solid red line). 




Figure 15. Antisymmetric solutions Eq. (|66p for c = 1; m = (solid black line), m = 2.726324772 
(solid blue line), m = 3.5 (dotted green line), m = 5 (long-dashed black line), m = 5.5 (dashed 
black line), and m = 5.875920066 (solid red line). 

Thus, although Mathieu functions have been sufficient to characterize a part of the spec- 
trum of the Schrodinger analog of our problem, we actually need to consider confluent Heun 
functions to cover all the cases. In other words, even when we achieved fully analytical 
solutions of the quantum analog differential equation, the spectra appeared just discrete not 
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Figure 16. Symmetric Heun solutions, Eq. (j65p . for m = 4.30206964 and several values of c: 
c = 0.7 (solid black line), c = 1 (long-dashed blue line), c = 1.5 (solid red line), c = 2 (dashed 
green line) , and c = 3 (dotted black line) . 




Figure 17. Antisymmetric Heun solutions, Eq. ()66p . for m = 2.72632477 and several values of c: 
c = 0.7 (solid black line), c = 1 (long-dashed blue line), c = 1.5 (solid red line), c = 2 (dashed 
green line), and c = 3 (dotted black line). 



revealing that some of them could be eventually continua. 



The set of confluent Heun functions therefore provide all the possible physical solutions 
of the actual problem in the 5D space. This was not apparent from the Hamiltonian point 



of view which assumes the Sturm-Liouville operator "H 
physical situation. 



to represent the 
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VI. FINAL REMARKS AND CONCLUSION 



In order to physically assess massive modes, one can evaluate the variation of the effective 
gauge coupling as a function of the Kaluza-Klein masses. Actually, KK contributions can 
not be significant as compared with the Coulomb potential because the coupling of massive 
modes to (fermion) matter on the brane develops a Yukawa type potential in the non- 
relativistic limit. To show that this is a decreasing function of m we should evaluate the 
coefficients that multiply the relevant sector of the four-dimensional action 

-Jdy e^(^)-^n(.)/2 |^^2^^^(^) + ^ ^2^^^y^^ I d^^f^j,^_ (67) 

However, in order to simplify this computation we can assume that the coupling with the 
brane takes place exactly on the 4D ordinary space-time, namely at y = 0. It is precisely at 
this value of y where the relevant physical effects should be much stronger. For simplicity 
let us consider the series of the quantum analog eigenvalues which serves as a discrete 
representative of the continuum. Thus, the effective 4D electrostatic potential would read 

7 + E^<(0)j (68) 

where gi, q2 are two test charges separated a distance r in ordinary 3D space and the Kaluza- 
Klein masses m are numbered with n in ascending order. See Fig. [T41 where the first u"ewen{y) 
modes are fully displayed, and Fig. ITSl where the first and the tenth modes are compared. 
See Fig. [19] to appreciate the first 10 values at the origin. This, together with the negative 
exponential factor, essentially decouples the massive modes from the physics on the domain 
wall. Far from the membrane, all massive modes become constants like the zero-mode is, 
and as a consequence the 5D phenomenology results completely modified from ordinary 4D 
electromagnetism. See e.g. Refs. for the study of this issue in the case of gravity. 



In this paper we have studied bulk and four- dimensional gauge propagation modes in a 
warped extra-dimensional space with a dilaton field. We have set up a sine-Gordon thick 
membrane which bounces at the extra-coordinate origin. A five- dimensional metric was 
dynamically generated consistently with the soliton brane and the dilaton background. In 
such a framework we studied the solutions of a five-dimensional gauge field. 
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Figure 19. Sequence of the first KK values of lieven(O) for c = 3/8 displaying the relative weights 
of the KK modes on the brane. 

First, we have found the exact quantum- mechanical analog of our original five-dimensional 
stringy problem. We have shown that the corresponding Schrodinger potential function is 
a quotient of simple second- and fourth-order polynomials that we could solve analytically. 
We next obtained the exact quantum-mechanical analog eigenspectrum and used it as a 
guide to analyze eventually the general solution. A localized zero-mode corresponding to 
the ordinary photon was guaranteed for a dilaton coupling constant above Aq. In general, we 
have found that the gauge-field dynamics are analytically given by confluent Heun functions 
which we have displayed for several representative cases. Furthermore, in contrast to the 
quantum analog results, in the general approach the mass of the gauge-field modes can be 
arbitrary for A G (Aq, Ai). In any case, we have shown that the Kaluza-Klein gauge spectrum 
is strongly attenuated on the brane as compared to the zero-mode of the theory. On the 
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other hand, we observed that in the bulk, far from the brane, the amphtude of an infinite 
tower of massive modes gets progressively relevant. Interestingly, the quantum-mechanical 
discrete mass eigenfunctions are completely decoupled in that region. 
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